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Pour en savoir plus Doc. J 1 065 


armi les nombreux problèmes de génie des procédés que rencontrent 
l'ingénieur et le pharmacien travaillant dans les industries chimiques, pétro- 
lières, pharmaceutiques, cosmétiques et agroalimentaires, l'écoulement d'un ou 
de plusieurs fluides à travers un milieu poreux fixe ou mobile tient une place 
prépondérante. Il suffit de citer les principaux procédés unitaires du génie des 
procédés (séchage, fluidisation, sédimentation, cristallisation, distillation, 
échange d'ions, extraction liquide-liquide...) pour voir le nombre pléthorique de 
canalisations, de colonnes, de cuves et de réacteurs au sein desquels le ou les 
procédés sont réalisés. 

De même, pour les industries de la santé, la formulation nécessite la concep- 
tion, le développement, la production et l'écoulement de matériaux poreux (ou 
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non poreux) fonctionnant par leur composition, leur préparation et leur agence- 
ment pour délivrer une action et rendre un service (par exemple la galénique). 

Après un bref rappel des principes fondamentaux de la mécanique des fluides 
appliqués aux cas d'écoulements de fluides parfaits ou visqueux newtoniens 
dans les conduites, ce texte fournit les notions de base indispensables sur 
l'hydrodynamique des écoulements dans les milieux poreux rencontrés dans les 

procédés de sédimentation et granulation, de réactions nécessitant un garnis- 
sage, de fluidisation et de filtration. 

Notre but est de proposer la ou les relations qui existent : 

— entre le débit de fluide et les propriétés caractéristiques du milieu poreux 
mobile et des fluides pour maintenir ce milieu poreux dans les conditions opti- 
males de fonctionnement afin de réaliser le procédé et d'élaborer le produit 
voulu (sédimentation, fluidisation) ; 

ou 

— entre les pertes de charge nécessaires pour assurer un débit connu et opti- 
mum de fluide, compte tenu des propriétés caractéristiques du milieu poreux 
fixe (réacteurs, filtration). 

Il est bien entendu que ce texte ne se veut nullement exhaustif et le lecteur se 
reportera utilement aux ouvrages hautement spécialisés présentés dans la 
bibliographie, pour une connaissance plus approfondie sur tel ou tel procédé. 


1.1 Fluides parfaits incompressibles 
en écoulement permanent 


On considère le fluide de masse volumique p initialement entre 
les sections A-|Bi et A 2 B 2 puis entre les sections A^B V A 2 A 2 au 
temps df. On désigne par l'aire de la section normale à l'axe A^ 
et par tq la vitesse dans cette section, A 2 l'aire en A 2 B 2 et par u 2 la 
vitesse dans cette section (figure 1). 


Figure 1 - Écoulement d'un fluide parfait incompressible 
dans une canalisation 



Les équations de bilan sont les suivantes : 

— conservation de la matière 

p A : u : =p A 2 u 2 = q m ( 1 ) 

avec q m débit massique du fluide (kg ■ s -1 ) ; 

— conservation de la quantité de mouvement (équation 
d'Euler) 


q m (u 2 -~up = ~F B 


( 2 ) 


avec F g résultante des forces extérieures (pesanteur, pression) ; 

— conservation de l'énergie (équation de Bernouilli) 


avec 


u? P, u? P 2 

_L + _L + z = JL + -1 + z ■ = c te 

2g pg 2 g pg 


g accélération due au champ de pesanteur (m 
P- 1, P 2 pression statique (Pa), 

Zi, z 2 cote, altitude (m). 



1.2 Fluides réels. Viscosité 


1.2.1 Définition 


Les seules forces qui existent dans un fluide parfait sont normales 
à la paroi. Dans les fluides réels que considère la mécanique des 
fluides interviennent en plus des forces de frottement ou de visco- 
sité, qui sont dans le plan de la paroi (fluides newtoniens). 

La vitesse u x des filets de fluide, parallèlement à l'axe Ox de la 
canalisation, varie suivant l'axe perpendiculaire Oz de la quantité 

du x 

par unité de longueur. 

dz 

Ce gradient de vitesse est accompagné d'une force s'exerçant 
entre deux filets voisins et dirigée suivant Ox. Cette force par unité 
de surface séparant deux filets voisins, appelée contrainte de frotte- 
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ment r, est proportionnelle au gradient de vitesse (relation de 
Newton) : 



Le coefficient de proportionnalité 77 , parfois noté /r, est la viscosité 

dynamique du fluide. On considère souvent aussi le quotient v = - 

P 

qui est la viscosité cinématique du fluide. 

Le tableau 1 rassemble les unités, dans différents systèmes, de 
ces deux grandeurs. 


Tableau 1 - Unités courantes des viscosités dynamique 
et cinématique 


Viscosité 

Dimensions 

SI 

CGS 

MKS 

V 



Poise 

Poiseuille 

(dynamique) 

M.L - 1 - T - 1 

Pa • s 

(Po) 

(dap) (1) 




10" 1 Pa ■ s 

10 Po 

v=n 



Stokes 

Myriastokes 

L 2 ■ T" 1 

m 2 ■ s" 1 

(St) 

(maSt) 

P 

(cinématique) 

1(T 4 m 2 - s -1 

10 4 St= 1 m 2 ■ s -1 


(1) dap : décapoise. 


1.2.2 Effet de la viscosité sur l'écoulement 


Dans une canalisation de diamètre d = 2 fî, on considère le fluide 
newtonien contenu dans un cylindre de rayon ret de longueur dx; 


p est la pression qui s'applique sur la face gauche et p+— dx celle 

dx 

qui s'applique sur la face droite (figure 2 ). 


En régime permanent, la vitesse u ne dépend que de la distance r 
à l'axe de la conduite. 

Un bilan des forces exercées sur le petit élément de volume rtr 2 dx 
s'écrit : 


1.2.3 Écoulement laminaire (f?e < 2000) 


Dans le cas où la vitesse d'écoulement est nulle lorsque le cylin- 
dre de fluide considéré a le même rayon que la canalisation (u = 0 
pour r- R), la vitesse u s'exprime : 


7- (R 2 ~r 2 ) ~ 
4?) 


dp'; 
d x) 


Le profil des vitesses est parabolique. 
Le débit volumique q v est : 


% 


pd/2 

Jo 


„ , 7Cd 4 ( 

u ■ 2 k- r ar = - 

1 2877V 


d P) 

dx) 


La vitesse moyenne, ou vitesse débitante, est : 


U = = _2Ü(_ dp' I 

n çfi 3277 V d x) 

4 

La chute de pression est donnée par la formule fondamentale de 
Poiseuille : 


dp 

dx 



et la perte de charge linéaire : 


(4) 


d H 32 77 
, — ~ 
dx pgd 2 


(5) 


Nota : Le signe « - » rappelle que les pressions décroissent au fur et à mesure que l'on 
progresse dans le sens de l'écoulement. 

Avec les termes adimensionnels : 


facteur de frottement - = 

2 4 14 dx 

P u m d 

et critère de Reynolds fie = 

V 

l'expression de la perte de charge (5) devient : 


f = _8_ 
2 fie 


Force de pression + Force de frottement = 0 


1.2.4 Écoulement turbulent (f?e > 3000) 


- Jtr 2 — dx- 2nr dx 77 — = 0 
dx dr 

d'où la relation fondamentale pour les vitesses : 

d u _ J_(_ dfn 

dr 277 v dxJ 



Figure 2 - Cylindre élémentaire de fluide newtonien en écoulement 


Profil des vitesses (profil de von Karman) 

On définit une vitesse fictive u* = u m 

• Au voisinage de la paroi, l'écoulement reste laminaire sur une 

U* V 

épaisseur y telle que < 5 . Dans ce cas : 

v 



u_ _ u* y 
u* v 


• Au centre de la canalisation, il existe un noyau turbulent qui 

U * y 

commence à une distance y de la paroi telle que > 30 . Dans ce 

v 

cas : 


— = 5,75 log ^+5,5 (7) 

u* v 
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( 7 * y 

• Dans la zone intermédiaire pour laquelle 5 < — - <30 , on a : 

v 

— = 11,5 log ^-3,05 (8) 

u* v 

Ces relations sont valables pour des conduites dites lisses. 

Pour des conduites rugueuses, où la hauteur moyenne e des 

/ u* e A 

aspérités dépasse l'épaisseur de la couche laminaire ^ > 5J , on 

remplace v * par e dans les équations (6), (7) et (8) donnant les vites- 
ses. 

Le tableau 2 présente quelques valeurs moyennes de hauteurs 
des aspérités e de canalisations industrielles en différents maté- 
riaux. 


Tableau 2 - Hauteurs moyennes des aspérités 
de canalisations industrielles 


Type de canalisation 

e (mm) 

Tube étiré (cuivre, plomb, inox, verre) 

0,001 à 0,010 

Tube en acier neuf 

0,01 à 0,10 

Tube galvanisé 

0,05 à 0,20 

Tube asphalté 

0,01 à 0,20 

Tube en fonte neuf 

0,8 à 1,5 

Tube en fonte oxydé 

0,8 à 1,5 

Tube en fonte incrusté 

1 à 3 

Ciment poli, bois raboté 

0,3 à 0,7 

Ciment brut, planches rugueuses 

0,7 à 2,5 

Maçonnerie 

5à 15 


Perte de charge 


On utilise l'abaque de la figure 3 représentant les variations de 


en fonction de Re pour différentes valeurs de rugosité relative - de 

d 

la paroi. 


Exemple 1 : soit une canalisation en acier galvanisé de 2" (tube 
gaz, série légère, NFE 29027) de diamètre intérieur réel d = 53,8 mm, 

e , 

de rugosité moyenne e = 0,15 mm (^ = 2,8 ■ 10 ). Cette conduite 

véhicule de l'eau à 20 °C (p = 1 000 kg • mr 3 ; p = 10~ 3 décapoi- 
ses). On considère successivement trois débits volumiques. 

• q v = 2,4 L min- 1 soit 4 • 10 -5 m 3 ■ s -1 . 


On peut d'ailleurs noter que, d'une façon générale, les écoulements 
laminaires conduisent à des pertes de charge très faibles. 

• q v = 16,5 L ■ min -1 soit 2,75 ■ 10 -5 m 3 • s -1 . 
u m = 0,121 m • s- 1 et Re = 6,51 ■ 10 3 

De plus, ^ = 2,8 ■ 10- 3 . Sur la figure 3, le point de rencontre de la 

verticale Re = 6,51 • 10 3 avec la courbe - intrapolée donne 

2 

- = 4,7 • 10- 3 . 

2 

On a donc f- — 1 = — — f/2 soit 5,2 ■ 10~ 4 m ■ rrr 1 . 

V dx J gd 

La perte de charge est de l'ordre de 5 cm de colonne d'eau pour 100 
m de conduite. 

• q v = 210 L min- 1 soit 3,5-1 0 -3 m 3 • s -1 . 
u m = 1,54 m ■ s- 1 et Re = 8,28- 10 4 . 

L'abaque de la figure 3 donne - = 3.4-10- 3 et il vient 

2 

(- — 1 = 6 ■ 1 0- 2 m ■ m- 1 , soit 6 m de colonne d'eau pour 100 m de 
V dx J 

conduite. 

Ces résultats sont résumés dans le tableau 3. 


Tableau 3 - Pertes de charge dans une canalisation 
en fonction du débit volumique (exemple 1) 


dv 

(m 3 ■ s -1 ) 

(m • s- 1 ) 

Re 

f 

2 

d H 
dx 

(mm 
d'eau ■ 
m- 1 ) 

(-d H) 
(mm d'eau 
pour 100 m 
de 

canalisation) 

4 ■ 10- 5 

1,76 ■ 10~ 2 

947 

0,00845 

CM 

1 

O 

CM 

2 

2,75 ■ 10~ 5 

0,121 

6 510 

0,0047 

52 ■ 10~ 2 

52 

3,5 ■ 10~ 3 

1,54 

82 800 

0,0034 

60 

6 000 


2.1 Mouvement d'une particule solide 
dans un fluide immobile 


On a u m = — — — = 1 ,76 • 10- 2 m • s -1 , d'où Re 
tid 2 /4 


Umdp 

1 


947 . 


Le régime d'écoulement est donc laminaire (Re< 2 ■ 10 3 ) . 


f 8 

On en tire _ = — = 8,45 • 10- 3 , valeur que l'on peut bien sûr aussi 
2 Re 

lire directement sur l’abaque de la figure 3. 

La perte de charge linéaire sera 


d H 
dx 


_ 4u 2 m f 
gd 2 


= 1,98- 1 0- 5 m ■ m- 1 . 


La perte de charge est donc de l'ordre de 2 mm de colonne d'eau 
pour 100 m de conduite. 


On considère une particule sphérique, de masse volumique p g , de 
diamètre d g , ayant une vitesse de sédimentation U, vitesse relative 
par rapport au fluide de masse volumique pf et de viscosité p dans 
lequel elle sédimente. 

On définit un critère de Reynolds de grain : 


fle„ = 


Pf Ud g 

n 


Le grain est freiné dans son mouvement par le fluide et il apparaît 
une force de traînée F qui dépend du critère de Reynolds et de 
l'aire A du maître-couple : 


F = Ne A 


p, U 2 
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Facteur de frottement 
J_ _ gd d H 
2 4 u m 2 dx 


Rugosité 
relative — 



Figure 3 - Variation du facteur de frottement 
f 

2 en fonction du critère de Reynolds Re pour 


Perte de charge dans une canalisation 


différentes valeurs de rugosité relative ^ 


avec Ne critère de Newton, appelée aussi C D , coefficient de traînée On utilise alors les critères adimensionnels X, proportionnel à d 3 , 

d'une particule. et V, proportionnel à L/ 3 : 


Le mouvement du grain obéit à des lois différentes suivant les 
valeurs du critère de Reynolds conduisant à différents régimes de 
sédimentation. Les expressions du critère de Newton Ne et de la 
vitesse de sédimentation U pour les différents régimes sont regrou- 
pées dans le tableau 4. 

Dans la quasi-totalité des cas, on connaît les masses volumiques 
du fluide et du grain, ainsi que la viscosité du fluide. Le problème 
posé est alors : 

— soit de calculer la vitesse de chute L/d'un grain de diamètre d g 
donné ; 

— soit de déterminer le diamètre d'un grain, ayant mesuré sa 
vitesse de sédimentation. 




y = 


Re„ 


3 JJ 2 
3 P 2 a U 3 


= K' U 3 


Ne 4tj ( Pg-Pf ) g 

Les valeurs limites de Xet Yen fonction du régime sont reportées 
dans le tableau 5. 

Donc, si d g est connu et U inconnu, on calculera X que l'on 
comparera aux valeurs limites 24 et 440 000. Si, au contraire, U est 

1 

donné et cL inconnu, on calculera Y que l'on comparera à — et 


24 


2 270. 


Tableau 4 - Expression de la vitesse de sédimentation en fonction du régime 


Régime 

Re g 

Ne 

Vitesse de sédimentation 


Stokes 

Re g< 1 

24 

Ne = — — 

Re g 

U _ d l (Pg-Pf) 9 

18 T] 

(9) 

Intermédiaire 

(Allen) 

1 < fie g < 10 3 

Ne = 18 ' 5 

Reo.e 

d g U «Pg-Pf) S>]°' 714 

6,55 77°-428 p0,286 

(10) 

Newton 

10 3 < Re g < 4 ■ 10 5 

Ne = 0,44 ± 0,02 

u pH (Pg-Pf) ai 1/2 

(11) 


Pf 
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Tableau 5 - Valeurs limites des critères adimensionnels 
X et Ven fonction du régime 


Régime X 

Stokes X < 24 


Y 



Intermédiaire 


Newton 


24 < X < 440 000 

440 000 


— < Y <2 270 
24 

Y> 2 270 



(a) à fond fermé (b) à fond ouvert 


Figure 4 - Sédimentation collective d'une suspension de particules 
dans une colonne 


2.2 Sédimentation d'une suspension 
de particules en régime de Stokes 


Tant que la teneur volumique p de la suspension en matières soli- 
des reste inférieure à 0,4, on considère que la suspension se com- 
porte comme un liquide newtonien dont la viscosité apparente est : 

774 = 77 ■ lO 1 ' 82 ^ 

et la masse volumique apparente : 

P A ~ P Pg + Il ~P) Pf 

Pour déterminer la vitesse de sédimentation U, on utilise la rela- 
tion (9) pour une particule, présentée dans le tableau 5 en rempla- 
çant p f par p A et 77 par r\ A . 

Par ailleurs, lorsque la sédimentation a lieu de façon discontinue 
dans un récipient fermé à sa base (figure 4 a), l'accumulation de 
grains dans le fond du récipient chasse un volume égal de liquide 
qui remonte avec une vitesse ascensionnelle u par rapport à la 
paroi. 

Par suite, la vitesse U, qui est la vitesse de la particule par rapport 
au fluide, sera égale à la vitesse W (vitesse de la particule par rap- 
port à la paroi fixe du récipient) augmentée de la vitesse ascension- 
nelle u : 


U = W + u 


L'égalité des volumes solides descendants et des volumes liqui- 
des ascendants se traduit par : 

P W= (1 -p) u 

En combinant ces deux équations, on a : 


W=( 1 -P) U 

Enfin, pour tout le domaine de variation de p, les résultats expéri- 
mentaux peuvent être représentés par le rapport ^ , fonction de 

d„ 

Re a , — et p (relation empirique de Richardson et Zaki) : 

d R 


w=o-py 


L'exposant r dépend du régime d'écoulement et du diamètre du 
récipient d R . 

Dans le cas limite d'un récipient de très grand diamètre (d R »d a ), 
on a : 

r = 4,6 en régime de Stokes, avec Re g « 0,2 
r= 2,4 en régime de Newton, avec Re g > 500. 


La variation de r en fonction de Re„ et de — entre ces deux limi- 

9 d R 


tes est donnée dans le tableau 6 . 


d 9 

Tableau 6 - Valeurs de r en fonction de Re g et de g- 

Re g 

Exposant r 

0 < Re g < 0,2 

r= 4,6 + 20 ^ 

d R 

0,2 < Re g < 1 

r= ( 4 , 4 + 18 —]Re-°' 03 
l d R ) s 

1 < Re g < 200 

r = [ 4,4 + 18 — Ifîe^ 0 ' 1 

l d R J 9 

200 < Re g < 500 

'= 4,4-fiegO. 1 

Re g > 500 

'st 

cm" 

II 


La comparaison théorie-expérience montre que l'accord est excel- 
lent uniquement en régime de Stokes, régime le plus courant à 
l'échelle industrielle. 

Quand la sédimentation est à fond ouvert ( u = 0), les particules 
sédimentent avec la vitesse U = W (figure 4 b). 

Pour dimensionner la hauteur Z de la tour de sédimentation utili- 
sée pour un procédé nécessitant un temps de séjour t s (échange 
d'ions, agglomération, enrobage, séchage...), il vient : 

Z=t,W 


J 1 065 6 


Toute reproduction sans autorisation du Centre français d'exploitation du droit de copie est strictement interdite. 
© Techniques de l'Ingénieur, traité Génie des procédés 



ELEMENTS DE MECANIQUE DES FLUIDES 


Un globule fluide se distingue d'un globule sphérique rigide par 
sa viscosité tj g qui n'est pas infinie et par la tension interfaciale glo- 
bule/phase continue a. 

Pour un globule sphérique de diamètre d, qui n'est soumis ni à 
des déformations ni à des oscillations, et en écoulement laminaire 
(Re g < 1), la vitesse du globule est égale à la vitesse de chute libre 
donnée par la formule de Stokes (9), multipliée par un facteur cor- 
rectif, coefficient de Hadamard H : 


H = 


3t? + 3 1 ] g 
2 t] + 3t] 9 


( 12 ) 


Les globules de grande dimension ne restent pas sphériques. 
Sous l'influence des forces normales de pression dynamique, ils se 
déforment. Cette déformation résulte de l'équilibre des forces dyna- 
miques, hydrostatiques et de tension interfaciale. Elle consiste en 
un aplatissement, ce qui augmente l'aire du maître-couple, donc de 
la traînée et par suite, une diminution de la vitesse de déplacement. 
De plus, une turbulence apparaît à l'arrière du globule. 

Comme en général les forces de viscosité interviennent en même 
temps que les forces interfaciales, on introduit trois critères adimen- 
sionnels Eôtvos, Weber et P, qui permettent de comparer les forces 
s'appliquant sur les particules : 


c .. 9 ( Pg-Pf ) d 2 , . . .. pesanteur 

Eo = , c est-a-dire — — — — p- 

(j tension interfaciale 


We 


P, u2d 


c'est-à-dire 


inertie 

tension interfaciale 


t? 4 9 (Pg-Pf) 


3.1 Vitesse de déplacement d'une goutte 


L'allure des variations de la vitesse d'une goutte en fonction de 
son diamètre, représentée sur la figure 5, fait apparaître l'existence 
de trois régimes, caractérisés par le critère de Froude Fr( rapport des 
forces cinétiques aux forces volumiques). 

Régime de Stokes 

La goutte est une sphère rigide. 


Fr 


P f U 2 

9d (Pg-Pf) 



Régime intermédiaire 

La déformation de la goutte n'est plus négligeable. 

La vitesse passe par un maximum quand le diamètre croît. 


Fr= 0,116 (Re g ) V2 



Régime d'Eôtvos 

La déformation est grande et la vitesse de la goutte est constante. 
Fr = 2,3 ■ EcT 1/2 



gouttes sphères rigides 

Figure 5 - Allure des variations de la vitesse de déplacement 
d'une goutte en fonction de son diamètre 



bulles sphères rigides 

Figure 6 - Allure des variations de la vitesse de déplacement 
d'une bulle en fonction de son diamètre 


3.2 Vitesse de déplacement d'une bulle 

Quand une bulle monte dans un liquide, l'allure des variations de 
sa vitesse en fonction de son diamètre, représentée sur la figure 6 , 
laisse apparaître l'existence de cinq régimes. 

Régime de Stokes 

La bulle est une sphère rigide. 

Régime de circulation interne 

La vitesse de la bulle peut être au maximum augmentée de 50 % 
par rapport à celle de la loi de Stokes (9), selon le coefficient H (12). 

Régime intermédiaire 

La déformation de la bulle n'est pas négligeable. La vitesse passe 
par un maximum. 
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Figure 7 - Écoulement d'un fluide à travers un lit de particules 


| 4 m 



Figure 8 - Milieu poreux. Modèle de Kozeny 


Régime d'Eôtvos 

La déformation est importante. La vitesse est indépendante du 
diamètre. 


4.2 Modèle de Kozeny. Relation 
de Kozeny-Carman 


Régime de Taylor 

Ce régime est obtenu pour les très grosses bulles en forme de 
champignon. On a alors : 

U =* soit Fr = 0,5 

V 2 


4.1 Loi de Darcy. Définition 
de la perméabilité 


Considérons une conduite de section droite S contenant un lit 
poreux de hauteur Z à travers lequel coule un liquide de viscosité tj 
à un débit volumique q v (figure 7). Darcy a observé que la perte de 
charge AP à travers le lit est proportionnelle à la hauteur Z du lit, à la 

q 

vitesse en fût vide, ou vitesse moyenne, u m = — et à la viscosité ry : 


u 


m 


B AP 
V Z 


(13) 


B étant la perméabilité dont la dimension est le carré d'une lon- 
gueur. 

Cette relation est parfaitement vérifiée en régime laminaire, la 
dégradation d'énergie provenant uniquement des frottements vis- 
queux. Si la vitesse croît, l'expression de la perte de charge se com- 
plique par l'introduction d'un terme proportionnel au carré de la 
vitesse moyenne, preuve d'une dégradation d'énergie par turbu- 
lence. 

Pour tenir compte de ces considérations et déterminer la perte de 
charge en fonction de la vitesse u m et des propriétés caractéris- 
tiques du milieu poreux, on utilise le modèle de Kozeny. 


Le milieu poreux réel, extrêmement complexe, est remplacé par 
un modèle simplifié, décrit mathématiquement à l'aide d'un nombre 
réduit de paramètres indépendants. 

Les propriétés caractéristiques du milieu poreux sont la 
porosités (le volume de vide du milieu poreux est sSZ) et la surface 
spécifique d'un grain a g constituant ce milieu poreux : 

surface d'un grain _ surface de tous les grains 
volume d'un grain volume de tous les grains 

On assimile le milieu poreux de volume V = SZa un faisceau de N 
pores identiques, droits et indépendants, de longueur Z p et de dia- 
mètre dp, mais dont la section n'est pas circulaire (facteur de circu- 
larité 7 = 1 ± 0,15). Ce modèle est représenté figure 8. 

Z D 

On appelle tortuosité le rapport T = -ÿ . 

On suppose que le régime d'écoulement est laminaire dans les 
pores. La relation de Poiseuille (4) s'applique donc : 


u 


mp 


dp AP 
32yil z p 


En écrivant : 

— l'égalité du débit de fluide en fût vide et du débit de fluide dans 


I ^ m Zp 

les pores u mp = — j I 


— l'égalité des surfaces internes du milieu poreux et du modèle 
a g V C\ - e) = N k dp Z p ; 

— l'égalité des porosités du milieu poreux et du modèle 

nd% 

sV = N — - Z , 

4 p 

la relation de Poiseuille (4) devient la relation de Kozeny-Carman : 


A P 
Z 


— = h K T] U m a g 


2 (1 -£) 2 


(14) 


La constante de Kozeny h K vaut : 

— pour des empilages de grains isométriques avec des porosités 

n'excédant pas 0,7 à 0,8, = 4,5 ± 1,5 et en pratique 5 ; 

— pour des empilages de fibres à porosité plus élevée, h K = 9 
pour e = 0,94 et h K = 22 pour e = 0,99. 
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La relation de Kozeny-Carman (14) peut s'exprimer en fonction du 

diamètre d a = — , diamètre de la sphère possédant la même sur- 

a 9 

face spécifique que le grain : 

Af = 36 h K lïüi 1^ 

Z K dl e 3 


4.3 Régîmes laminaire et turbulent. 
Relation d'Ergun 


On constate pratiquement que la valeur critique du critère de Rey- 
nolds de pores Re p au-delà de laquelle on quitte le régime laminaire 
est de l'ordre de 6. De façon analogue, en définissant un critère de 
Reynolds de pores modifié Re p , on quitte le régime laminaire pour 

Ré p >'\ , 


avec Re p 


P u mpdp 


et Re n = 


P“n 


(1 - e) a tj 


Rep 

4 T 


On démontre alors que la perte de charge correspondante est de 
la forme : 


AP 

Z 


3„ (1 -£) 


P u n 


(15) 


où h B = 0,3 est la constante de Burke-Plummer. 

La perte de charge totale s'exprime alors par addition de deux 
quantités : celle due aux pertes de charge par frottement visqueux 
contre les parois, relation de Kozeny-Carman (14), et d'autre part 
celle due aux pertes par turbulence, relation de Burke- 
Plummer (15) ; c'est la relation d'Ergun : 


AP 

Z 


JJSg (l -£) 2 


P a a (1 - e) 


(16) 


qui s'écrit encore sous forme adimensionnelle : 


r 

2 



avec 


2 


£^ AP 

a g (1 -£)p u 2 , Z 


On considère un lit de particules de diamètre d g , disposées sur 
une plaque poreuse. Ce lit est traversé par un fluide ascendant dont 
on augmente progressivement la vitesse u m en fût vide (figure 9 a). 

Pour un certain débit, la perte de charge dans le lit, de hauteur ini- 
tiale Z, et de porosité £,-, est la suivante : 


AP = Z, 36 h I, — 
dl 


(1 -£,) 2 


Lorsqu'elle est égale au poids du lit par unité de surface (1 Pa = 
1 N ■ rrT 2 ), c'est-à-dire : 


AP =Z/(p g -Pf)( 1 - ejig, 


Z 


°go$^o°°o|Pc 

n OQ-P f?0 O O 0 0 o 0 = 


Plaque poreuse 



Î<" 


€n A P 



® 


AP = pg AH 




A fluidisation initiale : la force 

correspondant à la perte de charge 
doit vaincre non seulement le poids 
du lit (point D) mais encore les 
forces de tassement qui assurent 
une cohésion entre les particules. 

B fluidisation en régime stationnaire. 

C fluidisation maximale et début 
de transport pneumatique 


Figure 9 - Fluidisation 


le lit commence à se soulever. C'est le seuil de fluidisation 
(figure 9 b, points A et D), caractérisé par la vitesse minimale de 
fluidisation : 


u 


min 


dg ef ( Pg-Pf)9 

36 h K ri ( 1 - £,) 


Au-delà de cette vitesse (figure 9 b, segment BC), les grains sont 
individualisés et sujets à des mouvements variés dans toutes les 
directions. Mais, pour une vitesse donnée, leur position moyenne 
reste statistiquement fixe ; le lit a donc une hauteur Zconstante, qui 
croît quand le débit de fluide croît : 

Z _ 1 - e, 

Z, ~ 1 - £ 

La perte de charge reste constante pendant la fluidisation, la poro- 
sité du lit fluidisé augmente avec la hauteur Z et avec la vitesse u m : 

“m = 1 zh _eL 

u min ef 1 - £ 

Si on augmente encore le débit, on va atteindre la vitesse corres- 
pondant à celle qu'auraient les grains s'ils tombaient en chute libre 
dans le fluide immobile (figure 9 b, point C). On aura alors entraîne- 
ment des grains qui seront éliminés dans le sens courant. C'est la 
vitesse maximale de fluidisation (qui correspond à la vitesse de 
sédimentation des grains) : 


U max 


d l ( Pg~Pf)9 
18rj 


On a alors : 


Umax = 10 0 -£/) 
u min ef 

qui oscille généralement entre 50 et 100. 

Cette représentation du lit fluidisé est idéalisée. Elle correspond à 
une fluidisation homogène dans laquelle les grains sont uniformé- 
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ment dispersés dans tout le lit, la porosité locale étant constante 

d'un point à un autre du lit = 1 1 . 

Vp f J 

Mais la plupart des opérations industrielles sont réalisées en 

phase gazeuse f- » O . Dans ce cas, la fluidisation devient hété- 

\p, J 

rogène. Certaines parties du lit sont immobiles, tandis que d'autres 
sont traversées par de grosses bulles gazeuses de vitesse uB et 
contenant une faible proportion de grains en suspension. Ces bulles 
viennent éclater à la surface du lit : c'est le bullage ou bouillonnage, 
et même le pistonnage si ces bulles occupent la section droite du lit. 


Pour mieux connaître le comportement d'une poudre vis-à- 
vis de la fluidisation hétérogène, Geldart a proposé une classifi- 
cation des poudres qu'il répartit en quatre classes en fonction 
de leur masse volumique et de leur diamètre moyen (figure 10 ). 

Cette classification est reprise par Castel dans Manutention des produits en vrac, 
A 9 301, traité L'entreprise industrielle. 

Classe A : les poudres fusantes. Elles se caractérisent par : 

— une fluidisation particulaire qui se traduit par une forte 
expansion du lit avec la croissance de la vitesse de fluidisation et 
l'apparition du bullage quand la vitesse du gaz est nettement 
supérieure à la vitesse minimale de fluidisation ; 

— une désaération lente (3 à 6 mm ■ s" 1 ) quand la fluidisation 
est interrompue : ces poudres restent longtemps fluides ; 

— une tendance au fusage lors de la vidange d'une trémie, 
c'est-à-dire un écoulement de type fluide à travers l'obturateur. 

Classe B : les poudres sableuses. Elles ont en commun : 

— une fluidisation agrégative avec une faible expansion du lit 
tant que la vitesse minimale de fluidisation, à laquelle apparais- 
sent déjà les premières bulles de gaz n'est pas atteinte ; ces bul- 
les sont exemptes de solide et grossissent rapidement par 
coalescence au cours de leur ascension ; 

— une désaération quasi instantanée à l'arrêt de la fluidisa- 
tion. 

Classe C : les poudres fines cohésives. Les forces de surface 
interparticulaires y sont du même ordre de grandeur que les for- 
ces aérodynamiques de traînée exercées par le gaz sur les parti- 
cules. Ces forces de surface peuvent avoir des origines 
multiples : attractions moléculaires, électricité statique, eau 
adsorbée ou adhésivité du produit. Ces poudres sont difficiles à 
fluidiser et sont le siège de renardages (formation de canaux de 
gaz) dans les lits de grand diamètre ou de pistonnages (alter- 
nance de couches de gaz et de particules) dans les lits de petite 
section. Leur fluidisation, quand elle est possible, ne peut se 
faire qu'avec l'appoint d'énergie mécanique (agitation ou vibra- 
tion). 

Classe D : les poudres granuleuses. Les particules y sont de 
grande taille (d g > 1 mm) ou de masse volumique élevée ; elles 
se caractérisent par une vitesse de fluidisation élevée (lits gey- 
ser) et une fluidisation turbulente avec explosion de bulles en 
surface du lit. 


Exemple 2 : Séchage de granules en lit fluidisé. Mise en fluidi- 
sation du lit. Domaine de fluidisation. 

On désire sécher des granules sphériques de diamètre 500 gm 
essentiellement constitués de lactose de l'amidon et pour une très fai- 
ble part de principe actif. La masse volumique de la matière sèche est 
de 1 500 kg ■ m . La teneur en eau des granules est n s = 4 kg d'eau 
par kilogramme de matière sèche et des mesures préliminaires indi- 
quent qu'ils suivent la loi d'hygroscopicité standard. 

Le séchage est effectué par de l’air à 50 °C dont l'humidité relative 
est 0,15. La viscosité de l'air à 50 °C est 1,9 ■ 10 ~ 5 dap et sa masse 
volumique dans les conditions normales de température et pression 
est de 1 ,292 kg ■ m -3 . 


La porosité du lit au repos est e = 0,4 . 


Déterminons tout d'abord la vitesse limite séparant les écoule- 
ments laminaire et turbulent de l'air dans le lit de granulés et la 
perte de charge correspondante 

L'écoulement laminaire cesse quand le nombre de Reynolds de pore 
modifié dépasse 1 : 


donc pour 


Re p = 


P “m 

(1 - e) a g t] 


(1 - e) a n 

U = - — 

m p 

Dans le cas des sphères, la surface spécifique de la particule est 

3=4 = 1,2 - 10 4 m- 1 . 
y d 

u 9 

On en déduit 


u 


m 


0,6 x 6 x 1 ,9 ■ 1Q - 5 
5- 10- 4 x1,3x^| 


12,55 10 - 2 m • s - 1 


La perte de charge résultante par unité de hauteur de lit est alors 
donnée par la loi de Kozeny-Carman (14) avec = 5 et 


(1 -e ) 2 _ ( 0 , 6) 2 
e 3 (0,4 ) 3 


5,65 . 


Si Z est la hauteur de lit, il vient donc 
A P 

y = 5x1,9- 10- 5 x 12,55- 10 ~ 6 x (1 ,2 ■ 1 0 4 ) 2 x5,65 

= 9 720 Pa ■ m - 1 


Comparons maintenant de cette perte de charge avec le poids du 
lit par unité de hauteur et par unité de section droite. 

Si S est la section droite de la colonne et Z la hauteur du lit, le poids 
du lit sera P= SZ (1 - e) p s g 

avec p s masse volumique du matériau humide. 

P ,, -v 
— = (1 -e) p s g 

La masse volumique du matériau humide est la moyenne harmo- 
nique des masses volumiques de la matière sèche et de l'eau pondé- 
rée par les titres massiques de la matière sèche et de l'eau dans le 
granulé : 


1 _ 1/14 0,4/1 ,4 _ 0,713 0,286 

Ps~ P M P B ~ 1 500 1 000 


0,760 ■ 10 - 3 m 3 - kg - 1 


P s ~ 1 300 kg • m - 3 

Le poids du lit par unité de section droite et unité de hauteur devient 

alors : — = 1 300 x 0,6 x 9,81 = 7 650 Pa - m - 1 . 

SZ 


Cette grandeur est inférieure à la perte de charge calculée précé- 
demment : au début de la fluidisation du lit, l'écoulement de l'air 
est donc laminaire. 


Calculons la vitesse minimale de fluidisation, en début de 
séchage. 

D'après les résultats précédents (régime laminaire), on peut calculer 
la vitesse minimale de fluidisation en utilisant le modèle de Kozeny : la 
perte de charge doit équilibrer le poids du lit par unité de section droite. 


— = 7,680 • 10 3 Pa - m - 1 = — = h K 77 3 2 (1 e)2 

Z 9 e 3 

7,680 


SZ 

Il vient u. 


5x1,9x1,44x5,65 


= 1 Q - 1 m • s - 1 
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Nota : supposons des grains dix fois plus gros, soit 5 mm de diamètre, toutes choses éga- 
les par ailleurs. La surface spécifique serait donc divisée par 10, la perte de charge par 100, le 
poids du lit ne variant pas. Le début de la fluidisation se passerait donc en écoulement turbu- 
lent et la vitesse minimale de fluidisation devrait être calculée non plus à partir de la formule de 
Kozeny-Carman (14), mais à partir de celle d'Ergun (16). 

Calculons enfin la vitesse d'entraînement en fin de séchage. 

Admettons que le grain soit complètement sec, sans avoir subi de 
variation de volume. L'eau étant remplacée par de l'air dans les pores, 
la masse volumique du grain (et non de la matière sèche) a diminué. 
Négligeons la contribution de l'air, 1 500 fois plus léger que le matériau. 

I m 3 de matériau pesait à l'origine 1 300 kg. La fraction pondérale 
0,297 d'eau lui est totalement enlevée. Le volume de 1 m 3 de matériau 
pèse alors 1 300 x 0,713 = 930 kg. La masse volumique du granulé 
devient alors p s = 930 kg • m- 3 . 

Pour connaître le régime d'écoulement, on calcule le nombre X : 
v 4 P (P S -P)9C(| 41,09x930x9,81 x (5) 3 x (1 0- 4 ) 3 „ cnn 

3 r, 3 3 (1 ,9) 2 x (10- 5 ) 2 

On se trouve donc en régime intermédiaire d'Allen (10). 

On sait que dans ce régime X = fie 3 Ne et Ne = L ... C? 

y Re g 

II vient donc fie’' 4 = — = 248 . 

9 18,5 

pUd n 

On en déduit fie„ = , 

b p 

il vient u = 1 ,80 ms 4 

On peut aussi calculer la vitesse en régime intermédiaire de Allen 
(10) qui donne li = 1,89 ms 4 . 

Il en résulte que l'on peut fluidiser entre les deux vitesses 
limites 

0,1 m ■ s 1 < U< 1,90 m ■ s 1 


6 000 

Pg-Pf 
(kg ■ m" 3 ) 

1 000 


100 

10 100 1 000 d s (gm) 

Figure 10 - Classification des poudres 
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Le but de la filtration est la séparation d'une phase continue et 
d'une phase dispersée, initialement mélangées ; ce mélange est 
désigné sous le terme général de suspension (ou préfiltre). 


La séparation peut être effectuée de deux façons différentes. 

Filtration dans la masse 

La suspension s'écoule à travers une masse poreuse (lit de parti- 
cules consolidé ou non, papier, feutre...) à l'intérieur de laquelle les 
particules ou globules de la phase dispersée sont plus ou moins 
retenus, entraînant un colmatage progressif de la masse filtrante. La 
phase continue sortant de la masse poreuse, appelée filtrat pour les 
liquides, sera soit parfaitement pure, soit appauvrie en phase dis- 
persée. 



Figure 11 - Filtration des liquides 


Filtration sur support 

On amène la suspension au-dessus d'un support (toile, feutre, 
membrane, grille...) sur lequel les particules vont se déposer sous la 
forme d'un gâteau d'épaisseur croissante. Le filtrat est recueilli sous 
le support. 


6.1 Lois générales 

On considère l'expérience de filtration d'une suspension schéma- 
tisée sur la figure 1 1 . 

À l'instant t, l'épaisseur du gâteau est Z et le volume du filtrat 
recueilli est V. 

Une pression AP est appliquée au-dessus de la suspension. 

Le problème posé par la filtration consiste à déterminer les varia- 
tions du volume de filtrat V, de la perte de charge APet de la hauteur 
du gâteau Z en fonction du temps t, compte tenu des propriétés du 
préfiltre et du gâteau. 

Appelons q v le débit volumique instantané de filtrat et posons 

q v =^=nu, (17) 

df 

£2 étant la surface filtrante. Supposons par ailleurs que le gâteau est 
homogène et incompressible et que le régime d'écoulement à tra- 
vers le gâteau est laminaire. 

On observe que la perte de charge AP est proportionnelle au débit 
instantané et à la viscosité du fluide (loi de Darcy, § 4.1) : 

AP = R pu (18) 

avec fi somme des résistances du support R s et du 

gâteau en cours de formation fig (m -1 ). 

On définit alors une résistance spécifique du gâteau a telle que 

fi G = a M/Q (19) 

avec M = nZ(1-e)pQ masse de gâteau séché, 

Pq masse volumique du gâteau, 

donc M/n masse de matériau séché déposé par unité de surface 
filtrante. 

La comparaison avec la loi de Darcy (13) où intervient la perméa- 
bilité 6 du gâteau donne : 


1 

a = 

B (1 - e) p G 
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6.1.1 Relation entre l'épaisseur Z du gâteau 
et le volume de filtrat V 


On suppose que la masse M du gâteau est à chaque instant pro- 
portionnelle au volume V de filtrat recueilli : 

M = W V 


avec W masse de particules ou de gâteau séché déposé 

par unité de volume de filtrat. 


La filtration d'une masse S de suspension donne : 

— un gâteau séché de masse M = s S ; 

— un gâteau humide de masse M' = m s S ; 

— une masse de filtrat pf l/ = (1 - m s) S, 


s teneur en particules de la suspension, 

pf masse volumique du filtrat, 

. ,,, ..... /'masse de gâteau humide^ 

m rapport d humidité 1 5 1 

nZO-e)p G p f s 


masse de gâteau séché J ' 


d'où l/V = 

V 1 -ms 

avec la relation entre le rapport d'humidité m et la porosité e : 


m = 


1 + 


ePf 

(1 -e) p G 


6.1.2 Vitesse de filtration 

En utilisant l'expression de la perte de charge (18) avec celles de 
la résistance du gâteau Rq (19) et de la vitesse de filtration u (17), il 
vient : 


AP 


r R OWVY dV 

v n Jn dt 


( 20 ) 


Si la filtration est à pression constante (sous vide ou sous pres- 
sion), alors A P= C te et l'équation (20) devient : 

f _ ^s 7 ! aWrjV 
V P2AP 2£i 2 A P 


6.1.3 Résistance du support R s et résistance 

spécifique du gâteau a. Équation de Ruth 


On définit : R s = aWV 0 /£2 

avec V 0 volume fictif de filtrat qui correspondrait à la formation 
d'un gâteau de même résistance que la membrane. Le 
temps nécessaire à cette formation est f 0 . 

Il vient, en utilisant l'expression de la vitesse de filtration (17) et en 
intégrant à pression constante 


avec a 


aWr\ 

2 Pi 2 A P 


fo + f = a ( V 0 + V ) 2 



Figure 13 - Filtre-tambour pour filtration continue sous vide 

L'utilité de cette équation réside dans la comparaison de la résis- 
tance spécifique du gâteau et de la résistance du support pour 
savoir s'il est raisonnable de négliger cette dernière. 



Figure 12 - Variation des pressions en fonction du débit 
pour une filtration à pression et débit variables 


6.1.4 Filtration à pression et débit variables 


C'est le cas des filtres alimentés par une pompe centrifuge dont 
on suppose connue la courbe caractéristique pression-débit. On 
connaît également la courbe donnant la perte de charge P i - P 0 du 
support et des canalisations du filtre en fonction du débit. Ces deux 
courbes se coupent au point A qui donne le débit initial de filtrat 
(figure 12 a). 

À tout instant de la filtration, la différence des ordonnées mesure 
la chute de pression P 2 - P-] à travers le gâteau seul : 


77 faWy\dy 
Pi V Pi J dt 


soit encore 


v = 


n 2 


r/al/)/ dV/dt 


Connaissant — - — , le volume de filtrat V est donc obtenu graphi- 
tai/)/ 

quement lorsque le point figuratif décrit la courbe caractéristique 
P 2 = f (débit) de A à B. 

Le temps nécessaire à l'obtention de ce volume est obtenu gra- 
phiquement (figure 12 b) par : 

r v 

t= f 

J 0 dV/dt 

Nota : cette représentation de t/V en fonction de l/est souvent utilisée pour détermi- 
ner expérimentalement les résistances spécifiques du support et du gâteau. 


6.1.5 Filtration continue sous vide 


Soit 1 jf = 9/2n la fraction de la surface totale Q du tambour à tra- 
vers laquelle a lieu la filtration à un instant donné et N la vitesse de 
rotation du tambour (figure 13). 

Le débit volumique de filtrat q v est obtenu par l'équation : 


Rs*? aWri 

Qv + — 

Pi AP 2 A/12 2 A P 


Q 


2 

V 


Le débit volumique de suspension q v est obtenu par l'équation : 

R s qsp 0 , ccs 2 pl 

¥ = Qv + Qv 

P2APW v 2 N£i 2 WA P 


avec s 
Po 


Po = 


Pf P a 

Pf ■ S+PgO -S) 
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